Evidence for short-range antiferromagnetic fluctuations in Kondo insulating YbBi 2 
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The spin dynamics of mixed-valence YbBi2 has been studied by inelastic neutron scattering 
on a high-quality single crystal. In the Kondo-insulating regime realized at low temperature, the 
spectra exhibit a spin-gap structure with two sharp, dispersive, in-gap excitations at fiu) ~ 14.5 
and ~ 20 meV. The lower mode is shown to be associated with short-range correlations near the 
antiferromagnetic wave vector qo = i, i). Its properties are in overall agreement with those 
expected for a "spin exciton" branch in an indirect hybridization gap semiconductor. 



Heavy-fermion compounds exhibit a whole spectrum of 
unconventional low-temperature behaviors, basically re- 
flecting the existence of a very small energy scale, of the 
order of a few tens of kelvin, in the electron subsystem 
. This energy scale is a hallmark of strong electron cor- 
relations and, in the metallic case, is associated with the 
Kondo temperature below which the heavy-quasiparticle 
Fermi liquid forms. In insulating compounds, on 
the other hand — so-called "Kondo insulators" (KI) 
or "mixed- valence semiconductors" (MVSC), such as 
CeNiSn, Ce 3 Bi 4 Pt 3 , SmB 6 , YbB 12 , or UPtSn — it cor- 
responds to the opening of a very narrow, temperature- 
dependent, energy gap in the electron density of states 
2]. The physical origin of this insulating state is still 
incompletely understood. It has been argued that a 
number of aspects can be explained in terms of a one- 
electron band picture, with a "hybridization gap" form- 
ing at low temperature in the electronic density of states 
at the Fermi energy Q. However, there is growing evi- 
dence that strong electron-electron correlations are cen- 
tral to the emergence of the gap behavior, and that their 
effects cannot be reduced to a mere renormalization of 
quasiparticle states. The spin dynamics of these systems 
is also peculiar: in most examples studied to date, in- 
elastic neutron scattering (INS) spectra typically exhibit 
a spin-gap response (A s ~ 1 — 10 meV) at low temper- 
ature, which seems directly related to the KI state and 
disappears rapidly when a single-site fluctuation regime 
is recovered by heating 0. Information on Q depen- 
dences has remained rather scarce, and to a large extent 
inconclusive, either because of complex anisotropy effects 
as in CeNiSn, or because measurements were carried out 
only on polycrystal samples. YbBi2 is a promising candi- 
date for further investigations: it is an archetype KI com- 
pound || with a simple NaCl-type crystal structure (in- 
terpenetrating fee sublattices of Yb ions and Bi 2 cuboc- 
tahedra), and previous inelastic neutron scattering (INS) 
experiments on powder 0, 0, have indicated the pres- 
ence of two narrow magnetic excitations near the spin- 
gap edge. Early single-crystal measurements [t| were in- 
terpreted in terms of a single dispersive low-energy mode 



with high intensity along [111], but the form of the Q de- 
pendence was not clearly established. In this Letter, we 
report a detailed investigation of the low-energy spin dy- 
namics in YbBi2 showing that there indeed exist two dis- 
tinct excitations, and that the lower one has a maximum 
in its intensity, together with a minimum in its energy 
dispersion, at the wavevector q = (|, s, |). We ascribe 
this behavior to short-range antiferromagnetic (AFM) 
correlations and argue that it could correspond to the 
spin-exciton branch recently predicted to occur in the KI 
regime of a periodic Anderson model [Icj . 

The sample studied consisted of two high-quality single 
crystals (total volume w 0.4 cm 3 ) grown by the traveling- 
solvent floating-zone method using an image furnace with 
four xenon lamps |ll| . INS experiments were performed 
on the thermal-beam triple-axis spectrometer 2T at the 
LLB (Saclay). The sample was oriented with a (110) 
crystal axis normal to the scattering plane, and cooled 
to temperatures between 10 and 80 K inside a closed- 
cycle refrigerator. Neutron spectra were recorded at fixed 
final energy, Ef — 14.7 meV (PG 002 monochromator, 
analyzer, and filter on the scattered beam), yielding an 
energy resolution of « 2 meV at zero energy transfer. 

As noted in earlier powder studies 0, S @; signif- 
icant nuclear scattering from acoustic and lower optic 
phonon modes exists in the energy window of interest 
15 < fiLO < 25 meV. Special care must therefore be ex- 
erted to ensure that this signal is correctly subtracted 
out from the experimental spectra. To this end, we have 
compared data obtained at several equivalent Q vectors, 
using lattice-dynamics measurements and numerical sim- 
ulations as a guideline. The procedure is exemplified 
in Fig. ^ for Q = T + {\i\i\), where r denotes a recipro- 
cal lattice vector. Experimental data are plotted together 
with the calculated nuclear and magnetic components. 
The decomposition was made by assuming that the Q 
dependence of the magnetic intensity follows the Yb 3+ 
atomic form factor, and that the phonon intensity varies 
proportional to Q 2 , as expected from the model calcula- 
tion for this particular set of Q vectors. The important 
check, shown in the middle frame of Fig. HI is that the 
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FIG. 1: (Color online) Separation of nuclear and magnetic 
scattering (see text) in spectra measured at T — 11 K for three 
equivalent L points, Q = Q(l, 1, 1); hatched areas represent 
the phonon signal, also plotted as dashed line. 



spectrum measured at the intermediate vector (|, |, |) is 
correctly reproduced using the partial (phonon and mag- 
netic) contributions derived from the Q = (|, |, |) and 
(|,|,|) spectra (Fig. ^ upper and lower frames) [27) . 
The contributions corresponding to the zone boundary 
X point, q = (0, 0, 1), were similarly estimated from the 
(1, 1, 2) and (1, 1, 4) spectra, yielding a phonon signal 
comparable to that derived above. An average of the 
two was thereafter taken to represent the phonon back- 
ground for all measured zone boundary spectra. On the 
trajectory between (|, |, |) and the (1, 1, 1) zone center, 
the phonon intensity was assumed to scale with Q 2 [28| . 
Along (£,£,1), flat low-energy optic modes existin g at 
about 24 to 28 meV throughout the Brillouin zone |l2( 
hamper the separation of the magnetic contribution in 
this energy range. 

Typical magnetic excitation spectra obtained from the 
above treatment are displayed in Fig. |2 for three differ- 
ent Q vectors on the zone boundary, Q = Q = (§,§,§), 
(1.35, 1.35, 1.8), and (1, 1, 2), (see sketch of Brillouin 
zone in Fig. |3J). In the upper frame, one notes a clear 
spin-gap region (no detectable signal below w 10 meV) 
followed by two narrow peaks, Ml and M2, at 14.5 and 
19 meV. These energies coincide with those found in 
the previous powder measurements, confirming that the 
peaks arise from the same magnetic excitations. Moving 
away from the [111] direction, Ml is strongly suppressed, 
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FIG. 2: (Color online) Magnetic excitation spectra measured 
for three different scattering vectors at T = 11 and 80 K 



while shifting to higher energies. It becomes practically 
undetectable in the c* direction (bottom frame). In con- 
trast, the M2 peak remains visible for all three Q vectors, 
showing only a moderate increase in energy between the 
L and X points. The solid lines in Fig.Elrepresent fits of 
the data using a scattering function S mag constituted of 
three spectral components, in analogy with Refs. 
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The shapes of the normed spectral functions 
Pi(E,Q,T) were chosen Gaussian for Ml (nearly res- 
olution limited) and Lorentzian for M2 as well as for the 
broad "40-meV" peak. The position and width of the lat- 
ter component, whose maximum lies beyond the limit of 
our experimental energy window, were fixed (except for 
the form-factor Q dependence) at values obtained in the 
powder experiments |g. This is justified because the ex- 
citation results primarily from single-site processes |l3j . 

The same procedure was applied to spectra measured 
at a number of scattering vectors, both inside and on the 
boundary of the r = (1, 1, 1) Brillouin zone, in order to 
determine the Q-dcpendence of the lower two excitations 
Ml and M2. The most notable result is that the inten- 
sity of Ml has a pronounced maximum at the L point, 
which also corresponds to an absolute minimum in its 
dispersion curve. This is demonstrated in Fig. where 
the energy and integrated intensity of Ml are plotted as a 
function of the reduced q vector for trajectories through 
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FIG. 3: Q dependence of the energies and intensities of the 
low-energy peak Ml at T = 11 K; left frame: intensity map 
over one quadrant of the Brillouin zone (circles denote Q vec- 
tors at which energy spectra have been measured, with di- 
ameters indicative of intensities); right frame: intensities and 
dispersion curves along the zone boundary (circles) and the 
[111] direction (squares). 

the L point respectively parallel to [111] and along the 
zone boundary. The intensity drop when shifting away 
from Qo is accompanied by a moderate increase in the ex- 
perimental line width. An intensity map of Ml covering 
one quadrant of the projected Brillouin zone is presented 
in Fig. 01 K clearly shows the maximum centered around 
Qo, with a steep drop in intensity on going toward the 
zone center. Ml retains significant intensity all along 
possibly going through a local maximum near 
Q = In contrast, it is drastically suppressed 

along the entire (1, 1, C) direction. 

The corresponding analysis for M2 is less straightfor- 
ward because of the above-mentioned overlap with low- 
energy optic phonons. The general trend is a weaker en- 
ergy dispersion in comparison with Ml, and a much more 
uniform azimuthal distribution of the intensity. However, 
as in the case of Ml, a pronounced reduction of the in- 
tegrated signal occurs near zone center. The result is 
unambiguous because, in this region, M2 becomes nar- 
row enough and thus can be easily separated from the 
phonon background. In the numerical fits, it turns out 
that the loss in spectral weight results mainly from a 
decrease in the signal line width. 

With increasing temperature, one observes a com- 
plete suppression of the Ml excitation (see spectra for 
T = 80 K in Fig. in agreement with the powder re- 
sults |E @- The data also confirm the recovery of 
magnetic intensity, ascribed to quasi-elastic scattering, 
in the former spin-gap region. It should be noted that 
the disappearance of Ml takes place without apprecia- 
ble shift in the excitation energy. The T dependence of 
the parameters for Ml and M2 derived from the fit is 
summarized in Fig. 0] for Q = Qo- The Ml peak is 
steeply suppressed between 25 and 60 K, with an attend- 
ing increase in its line- width (left frame). On the other 
hand, M2 (right frame) is little affected by the increase 
in temperature up to 80 K. 
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FIG. 4: Temperature dependence of fit parameters for the 
low-energy peaks Ml (left) and M2 (right) in the energy spec- 
tra measured at Qo = (§,§,§); lines are guides to the eye. 



The main outcome of the present measurements is the 
unambiguous observation of two distinct excitations in 
the magnetic neutron spectrum at T — 10 K, exhibit- 
ing contrasted Q and T dependences. The lower one, 
Ml, is of particular interest to us because it is specific 
to the low-temperature KI state. The distribution of in- 
tensity over the Brillouin zone, showing a pronounced 
maximum near Qo = Tm + qo denotes the importance 
of magnetic fluctuations associated with the wave vec- 
tor qo = (■§, s, h). Short-range AFM correlations are 
central to the physics of metallic HF systems ^4|, es- 
pecially near quantum critical points, and there is evi- 
dence that they pl ay a role in KFs as well [l5L ITtij . It 
has been reported [Tj| that some metallic heavy rare- 
earth -RB12 compounds order in complex, amplitude- 
modulated, incommensurate structures. The paramag- 
netic ground state is specific to YbBi2, likely due to a 
Kondo spin- liquidlike phase [l|. The magnetic couplings 
responsible for the observed AFM correlations may be 
affected by Fermi surface changes due to the gap forma- 
tion. The steep decrease in the intensity of Ml taking 
place between 25 and 60 K is directly correlated with the 
filling of the gap observed over the same temperature 
range in optical conductivity spectra cr(u>) [18j. Basing 
on a parallel between the electronic structures of YbBi2 
and metallic MV YbAl3, Okamura et al. 0] have argued 
that the shoulder occurring below 40 meV in their low- 
temperature measurements of <r(ui) reflects the indirect 
gap in the renormalized density of states expected from 
a periodic Anderson model. In the neutron spectra, this 
energy of 40 meV roughly corresponds to the position of 
the broad peak in the time-of-flight data which disap- 
pears on heating in the same temperature region (T > 50 
K) as the anomaly in o~(uj), and we thus believe that it 
represents a plausible estimate of the spin gap value. If 
this is right, Ml has to be considered an in-gap mode, 
analogous to an exciton in a conventional semiconductor. 

This view is in line with the model proposed by Rise- 
borough in Ref. which predicts a sharp branch of 
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in-gap dispersive "spin-exciton" excitations to form, un- 
der certain conditions, in KI systems as a result of AFM 
exchange interactions J(Q). This process is favored by 
the enhancement of the real part of the dynamical sus- 
ceptibility x'(Q,cj) taking place at the spin-gap edge for 
Q vectors close to a zone boundary AFM point 29] . The 
excitations are described as the continuation of an AFM 
paramagnon branch which softens at low temperature 
and undergoes strong narrowing due to the suppression 
of the electron-hole decay channel when it eventually falls 
below the threshold of the Stoner continuum. Moving to- 
ward the zone center, J(Q) decreases and x'(Qi w ) broad- 
ens, leading to a positive dispersion, and the gradual dis- 
appearance, of the spin-exciton branch as observed in the 
present experiments. 

In this interpretation, the damping, and subsequent 
disappearance of Ml with increasing T is a direct conse- 
quence of the gap filling evidenced by optical and trans- 
port measurements. The model implicitly assumes that 
the gap forms through a hybridization process taking 
place coherently at each 4/ ion site, and would thus be 
expected to be sensitive to disorder caused by chemi- 
cal substitution on the Yb sublattice. However, neu- 
tron spectra for the Ybi_ x Lu x Bi2 solid solutions 
indicate that the spin gap is quite robust and can be 
traced up to 90% Lu dilution, which is at variance with 
its rapid suppression with increasing temperature (recov- 
ery of quasielastic scattering). This raises the interesting 
possibility that the spin gap could be driven by single- 
site effectsflocal spin-singlet state) as suggested by other 
authors [2ll Eij ] , without necessarily requiring the exis- 
tence of an insulating gap. Whether some degree of local 
character can be reconciled with the spin-exciton picture 
remains as a question for future investigations. 

In summary, we have demonstrated a significant ef- 
fect of AFM short-range correlations in YbBi2, leading 
to the formation of a sharp dispersive in-gap mode. The 
main properties of this excitation can be captured by 
a rather simple "spin-exciton" model, though the pos- 
sibility of local effects, as inferred from substitution ex- 
periments, should also not be overlooked. Low-energy 
peaks have so far been reported for a few KIs such as 
SmB 6 0, CeNiSn 0j, and now YbBi 2 , but failed to 
be observed in a number of others ■ It is noteworthy 
that mechanisms reminiscent of the spin-exciton forma- 
tion, but based on the existence of a superconducting 
rather than insulating gap, have been invoked to explain 
anomalous magnetic modes occurring in high-T c cuprates 
(the well-known "resonance peak") [25j or the HF com- 
pound UPd2Al3 j2||. The picture emerging here for an 
archetype KI may thus prove of interest for a broader 
class of "spin-gap systems" . 
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